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We study out-of-equilibrium dynamics of intense non-abelian gauge fields. Generalizing the well- 
known Nielsen-Olesen instabilities for constant initial color-magnetic fields, we investigate the impact 
of temporal modulations and fluctuations in the initial conditions. This leads to a remarkable co- 
existence of the original Nielsen-Olesen instability and the subdominant phenomenon of parametric 
resonance. Taking into account that the fields may be correlated only over a limited transverse 
size, we model characteristic aspects of the dynamics of color flux tubes relevant in the context of 
heavy-ion collisions. 



I. INTRODUCTION AND OVERVIEW 

Out of equilibrium dynamics of non-abelian gauge 
fields provides an outstanding challenge in quantum field 
theory. Its most prominent application concerns the de- 
scription of relativistic heavy ion collisions. It involves 
the solution of an initial value problem in real time 
based on the underlying theory of quantum chromody- 
namics (QCD). This is not accessible to standard lattice 
gauge theory simulations in Euclidean space-time. How- 
ever, important aspects of nonequilibrium gauge field dy- 
namics may be described using classical-statistical lat- 
tice gauge theory techniques in real time. Classical- 
statistical approximations to the quantum dynamics are 
expected to be valid if the expectation value of the anti- 
commutator of fields is much larger than the commuta- 
tor. Loosely speaking this is realized in the presence of 
sufficiently high occupation numbers per mode or large 
coherent field expectation values. 

The physics of nonequilibrium instabilities provides 
an important example where classical-statistical sim- 
ulations techniques can give accurate descriptions of 
the underlying quantum dynamics. This has been 
tested to high accuracy for self-interacting scalar quan- 
tum field theories [l|-l3| and Yukawa-like theories with 
fermions [4], where appropriate far-from-equilibrium re- 
summation techniques are available for the respective 
quantum field theory. In recent years classical simulation 
techniques have become an important building block for 
our understanding of the time evolution of instabilities 
in non-abelian gauge theories in the context of collision 
experiments with heavy nuclei 

The present work investigates classical dynamics of co- 
herent non-abelian gauge fields, which is also motivated 
by the notion of 'color flux tubes' that may form af- 
ter the collision of two Lorentz contracted heavy nuclei. 
These are characterized by intense color-magnetic as well 
as color-electric field configurations in longitudinal direc- 
tion. They are smooth in this direction and correlated 
over a transverse size associated to the inverse of the 
characteristic momentum scale Qs in the saturation sce- 



nario . To understand the time evolution of this con- 
figuration in QCD is a formidable task, which is further 
complicated by the longitudinal expansion of the system. 

In order to approach this complex question of non- 
linear gauge field dynamics, it is very instructive to con- 
sider first an extreme simplification which can also give 
analytical insights. For the ease of a later numerical 
treatment, we use SU{2) gauge theory. We neglect ex- 
pansion and employ a constant color-magnetic field con- 
figuration, with color index a = 1,2,3 and Lorentz 
index /i, whose spatial components i = x,y, z point in 
longitudinal (i = z) direction 



(1) 



Here the latter can always be arranged to point into the 
direction a = 1 in color space by a suitable gauge trans- 
formation without loss of generality. In the context of 
heavy ion collisions this configuration has been exten- 
sively discussed in Ref. [l3|, where the field in infinite 
volume is taken to be generated from the vector poten- 



tial Af. as 



Al^-\vB, Al^\xB, 



(2) 



with all other spatial components vanishing. Such a field 
configuration is known to exhibit a Nielsen-Olesen in- 
stability [ll| characterized by an exponential growth of 
fluctuations with maximum growth rate 



7NO = 



(3) 



where y/gB may be taken to be of order Qs- This expo- 
nential growth leads to a production of gluons, which is 
much faster than any conventional production process. 
Its consequences for the question of tlicrmalization in 
a heavy ion collision can be significant. Of course, in 
this context a single macroscopic constant color-magnetic 
background field is certainly not a realistic option and it 
is important to understand the robustness of the under- 
lying physical processes against suitable generalizations. 
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In a first generalization we will allow for temporal mod- 
ulations of the color magnetic background field configu- 
rations. This can be conveniently achieved by the gauge 
field configuration 

Al = 4 ^ (4) 

with all other components zero. It can be directly verified 
that this vector potential configuration gives the same 
form of the longitudinal magnetic field (HJ. However, 
there are important differences as compared to the pre- 
viously considered case, where the magnetic field is gener- 
ated from space-dependent vector potentials. In contrast 
to ([2]), the configuration (U) contributes to the non-linear 
part of the field strength tensor 

i^;, [A] = d^Al -d,A'l+ ge'^'^A^A: (5) 

from which the magnetic field is obtained as Bj = 
^tjkpa^^ Generating the magnetic field with 1^ via the 
non-linear term of the field strength tensor leads to co- 
herent oscillations in time with characteristic frequency 
~ y/gB. In contrast, it can be readily verified that ([2]) 
is a time-independent solution of the classical Yang-Mills 
equations. Accordingly, also the magnetic field derived 
from ^ is constant in both space and time as long as it 
is unperturbed. 

We discuss analytic solutions for the classical time evo- 
lution of the configuration (|4]) in temporal (Weyl) gauge 
and study the behavior of linear perturbations on top 
of the oscillating background field in Sec. HIl It is shown 
that these perturbations still exhibit robust growth sim- 
ilar to the Nielsen-Olesen-type. In order to go beyond 
the linear analysis we employ classical-statistical lattice 
gauge theory simulations in Sec. IIIIl The linear analysis 
is seen to accurately reproduce the numerical data for 
times which are short compared to the inverse character- 
istic growth rates. Afterwards, non-linear contributions 
become crucial. An important observation is that the 
non-linearities lead to a very efficient growth for all com- 
ponents of the gauge potential such that the details about 
the initial configuration ^ become irrelevant as long as 
sufficient primary growth is triggered. This adds to the 
robustness of the observed phenomena also for further 
generalizations of the initial configurations. 

An important step towards more realistic scenarios 
concerns the inclusion of fluctuations into the initial con- 
ditions. The initial configuration Q exhibits a sharply 
defined macroscopic field amplitude B which is certainly 
unrealistic, in particular, since no macroscopic colored 
objects can exist. Here the time-evolution of single 
macroscopic B-field configurations may be viewed as de- 
scribing the dynamics of individual members of an ensem- 
ble. Ensemble averages are then obtained from sampling 
initial conditions according to given averages or expec- 
tation values at initial time. In Sec. IIV Al we consider 
nonequilibrium ensembles, where we choose the initial 



values of the spatially homogeneous fields in Q randomly 
from a Gaussian distribution with zero mean. Again one 
finds robust growth also for the ensemble of homogeneous 
initial fields. A dramatic difference is that for the en- 
semble average one observes isotropization of the stress- 
energy tensor happening on a time scale much faster than 
l/\/gB because of a dephasing phenomenon. We com- 
ment on the possible relevance of this observation for the 
outstanding question of rapid isotropization in heavy ion 
collisions in Sec. IIVI 

Finally, we take into account that the fields should be 
correlated only over a limited transverse size, which in 
the color fiux tube picture is associated to the inverse of 
the characteristic momentum scale Qg. For this we divide 
in Sec. IIV 5] the transverse coordinate plane into domains 
of equal spatial size. Each domain is then filled with a 
coherent color magnetic field, however, using a different 
random color direction in each patch. If the transverse 
size is taken to be Qs or ^/gB then this should model 
characteristic aspects of the dynamics of color flux tubes. 
It turns out that one observes exponential growth for all 
components of this inhomogeneous gauge field configura- 
tion. The nature of the underlying instability still shows 
properties reminiscent of the Nielsen-Olesen-type. In 
particular, the characteristic dependence of the primary 
growth rate as a function of momentum still reaches its 
maximum at low momenta. In contrast, Weibel instabili- 
ties show a vanishing growth rate at zero momentum . 
This difference has also been emphasized in Ref. [l^l . We 
conclude in Sec. |Vl In an appendix we describe that the 
phenomenon of parametric resonance leads to a subdom- 
inant instability band at higher longitudinal momenta. 



II. LINEAR REGIME 

We first discuss analytic solutions for the classical time 
evolution starting from the coherent field configuration 
(ID). It is shown that linear perturbations on top of the 
oscillating background field still exhibit robust growth 
similar to the Nielsen-Olesen-type. We consider a non- 
abelian gauge theory with SU{2) color gauge group. Tak- 
ing into account two colors simplifies the analysis as com- 
pared to the SU{3) gauge group relevant for QCD and 
the difference is expected to be of minor relevance for 
the physics considered here [H, [13] . With the covariant 
derivative 

D';,''[A] ^ d^S^' + ge'^^'A';, (6) 

the classical equations of motion read 

iD^[A]F''''[A]r ^ 0. (7) 

In a first step, we will compute the classical time evo- 
lution starting from the coherent field configuration 
For an analytical understanding of the dynamics it is 
convenient to split the gauge field potentials A'^{x) into 
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a time-dependent background field A°(a-°) and a pertur- 
bation: 



Al,{x) = Al{x')+5A1{x) 



(8) 



We can then compute the time evolution in an expansion 
in powers of the perturbation (5v4^(a;). At zeroth order 
we obtain the field equation for the background field 



(9) 



The next order corresponds to the linearized equation for 
the fluctuations [T7| , 

{D^[A\D^'[A]5A''Y - {D^[A]D''[A]5A^^Y 
+ge''^''5AlF'""'[A] = 0. (10) 

Accordingly, equations ([9]) and (fTO|) correspond to the 
classical field equation ([7]) up to corrections of order 
{5Af. 

We choose temporal (Weyl) gauge where = 0. Writ- 
ing t = we consider the field configuration 



< 




t(Bg)' 



FIG. 1. Time evolution of the background field A{t) normal- 
ized to its initial value. It is d escribed by the Jacobi elliptic 
function cn{y/ Bg t, 1/2) \/B /g as a function of time t in units 
of the characteristic scale y/gB. 



Af{t) = A{t) {S-^S,, + S-H^y) 



(11) 



which corresponds to (U) at initial time with A{t = 0) = 
y/B/g for given initial color-magnetic field B. With the 
further initial condition dtA{t = 0) = the background 
field equation ([9]) reads 



We concentrate on momenta along the z-direction, 
i.e. we evaluate the modes of SAf{t, p) at Px ^ Py = 0. 
The evolution equation (fTO)) for the different compo- 
nents of SAf{t,pz) can then be written with SAf{t,Pz) = 
dtSAf{t,pz) as the following independent sets of coupled 
differential equations: 



d^A{t)+g'AitY = 



(12) 



and the solution is given in terms of a Jacobi elliptic 
function: 



Ait) 




(13) 



Here cni^/gB t, 1/2) is an oscillatory function such that 
A{t) = A{t + Ats) with period 



Ate = 



AK{l/2) 

VgB 



7.42 
TP' 



(14) 



g^A-^ + pi 
-igApz 



igApz 
g^A^ 



5AI 
5Al 



(16) 



The complex conjugate of the components {5A]^,5A'l) 
obey the same equation, which is obtained from (jlpj) with 
the replacement (Ml, M^) -^. {5Al*,6Al*). Similarly, 
one finds 



5AI 
5AI 



and 



2g2^2 



2g2A2 

2 42 , 



g'A 



^Al 
SAl 



(17) 



where if (1/2) denotes the complete elliptic integral of the 
first kind [1^ . Similar behavior has been discussed in the 
context of classical chaos |l5l| , or for parametric resonance 
in scalar field theories jl6|. The solution ([T^ is plotted 
in Fig. [T] For later use we note that the corresponding 
characteristic frequency is lub — 2ti/ /SiB — 0.847-^^5. 
We emphasize that no constant non-zero solution exists 
for A in this case. This is in contrast to the analysis of 
Ref. [l3], where the same gauge potential configuration 
in the presence of suitable external charges is assumed to 
lead to a time-independent A. 

In the following we insert this background field behav- 
ior into the linearized fluctuation equation (1101) and dis- 
cuss solutions for modes of 5A°i in spatial Fourier space: 



M?(i,p) = / d^xe-^P'^'6A1{t,^) 



(15) 



5AI 
5AI 
5A\ 




-g'^A^ -igApz 
pI_ igApz 
-igApz 2g^A^ 



SAl 
6AI 



SAl 



(18) 



Each of the linear differential equation matrices (|16l) - p^ 
with time dependent background field A{t) may be fur- 
ther analyzed by diagonalization. The time-dependent 
eigenvalues of the equation matrices read for (IK 



cot{pzr^g'A' + ?^±^x/49U^pl 



'Pi, 



for (HID 



i^2{Pzf 



■ig^A^ 



-P\ 



(19) 



(20) 
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pI), 
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and for ^ 

uJA{p.f = -[g^A^-pl), (21) 
^tiP^f = \ {^9^A^+pI ± ^Ig^A^ + Qg^A^pl+pt^ . 

The corresponding eigenvectors depend, in general, on 
A{t) and thus on time. However, the eigenvectors 
associated to u}2{pz)^, ^^iPz)"^ and ^4(^2)^ are time- 
independent and given by (l,l)/\/2, (— l,l)/\/2, and 
(l;l,0)/\/2, respectively. These include, in particular, 
the only negative eigenvalues 0^3(^2)^ — uJ4{pz)'^ for mo- 
menta p1 < g^A(t). These will play a particularly impor- 
tant role in the following since they turn out to govern 
the fastest time scales. 

It is very instructive to consider the example of ([T7| in 
diagonalized form, i.e. 

SA+=-{3g^A^+pl)SA+, (22) 
SA^^{g^A^-pl)SA^, 

where SA+ = 6Al + SA^ and SA^ = SAl - SAl. (The 
same equation for 5 A- would also be obtained from ^TE\\ 
by associating 6 A- to SA^+SAy.) These equations are of 
the Lame- type [l^: Since the squared background field 
appearing in ((22|) has periodicity Ais/2, each solution 
can be written as a linear combination of the form 

SA+(t + At/2,pz) = e^^+f"^) SA+{t,pz) (23) 

such that 

6A+{t,pz) = e2'^+(f^)*/^*- Tl+{t,pz) (24) 

with periodic functions Il+{t + At/2,pz) ~ Il+{t,pz) and 
similarly for 6A^{t,Pz)- The Floquet function C+{pz) 
is time-independent and leads to oscillating behavior if 
C+(pz) is real or to exponentially growing or decaying 
solutions if imaginary. 

The dominant exponentially growing solutions arise 
because of the appearance of the time-dependent neg- 
ative eigenvalues io^ipz)"^ = Lo^ipzY and we will concen- 
trate on them in the following. Subdominantly grow- 
ing modes are associated to parametric resonance and 
are discussed in detail in the appendix. We first ob- 
serve from the evolution equation (|22p for 5A^ {t, pz) that 
replacing the background field by a constant, A{t) -> 
\/B(t = 0)/g, would lead to the well-known Nielsen- 
Olesen result for the growth rate of modes with p^ < 
gB{t - 0): 

7No(p.) = ^gB{t = Q)-pl. (25) 

This is in accordance with the fact that the value for 
the maximum growth rate for constant fields, as stated 
in ([3]), is obtained for vanishing momenta. In contrast, 
since the background field is oscillatory in our case there 
will be deviations from the Nielson-Olesen result (1^51) . 
However, it turns out that to rather good accuracy their 



growth rates follow the Nielsen-Olesen estimate if the 
temporal average of the oscillating magnetic background 
field is used, which is explained in the appendix. The 
time average over one period AtB/2 of the square of the 
time-dependent background field ([T3|) . which enters the 
evolution equation for the fluctuation 6A^ (t, pz), is given 
for the above initial conditions by 

r(3/4)2 

= /}^/ , ^ gB(t = 0) 

r(5/4)r(i/4) ^ ^ ' 

M0A57gB{t = 0). (27) 

Indeed, replacing gB{t = 0) in (^5)) by the average value 
([27| reproduces the full numerical results in the linear 
regime at suSiciently early times to good accuracy, which 
is described in the following. 

III. CLASSICAL LATTICE GAUGE THEORY 

The exponentially growing solutions observed in the 
previous section will finally lead to non-linear behav- 
ior. In order to be able to describe this physics, we 
use classical-statistical Yang-Mills theory in Minkowski 
space-time following closely Ref. Q, to which we refer 
for further technical details. The above linear analysis 
will be seen to accurately reproduce the numerical data 
for times which are short compared to the inverse charac- 
teristic 'primary' growth rates of the linear theory. After- 
wards, non-linear contributions become crucial. In par- 
ticular, the primary growth induces 'secondary' growth 
rates which are multiples of the primary ones similar to 
what has been observed in Refs. Q, where no coherent 
magnetic fields were employed. As a consequence, a very 
efficient growth for all components of the gauge potential 
is observed such that the details about the initial configu- 
ration ([l} become irrelevant as long as sufficient primary 
growth is triggered. 

For the simulations we employ the Wilsonian lattice 
action for SU(2) gauge theory in Minkowski space-time: 

S[U] - -/3o E E { ^ U.,o^ + Tr C/I,o.) - ij 
E E { ^ (Tr U.,, + Tr - l| , 

X i<.j ^ 

(28) 

with X — (x", x) and spatial Lorentz indices i,j = 1, 2, 3. 
It is given in terms of the plaquette variable Ux.fiu = 

^l+o^^Pkvi where U^^^^ = Ux,^,^, ■ Here Ux,^, 
is the parallel transporter associated with the link from 
the neighboring lattice point a; -I- /I to the point x in the 
direction of the lattice axis = 0, 1, 2, 3. The definitions 
= 27Trl/5g and Ps = 2Trl/(g^-y) contain the lattice 
parameter 7 = as/ at, where denotes the spatial and 
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FIG. 2. Time evolution of spatial Fourier modes of as a 
function of time, all in appropriate units of the energy density 
e. The momenta are chosen in the z direction. 



at the temporal lattice spacings, and we will consider 
go = gs = g. The dynamics is solved in temporal axial 
gauge. 

Varying the action (f28|) w.r.t. the spatial link variables 
Ux,j yields the classical lattice equations of motion. Vari- 
ation w.r.t. to a temporal link gives the Gauss constraint. 
The coupling g can be scaled out of the classical equa- 
tions of motion and we set g — 1 for the simulations. We 
define the gauge fields as 

gAtix) = -^Tr{a^U,{x)) (29) 

where a^, and denote the three Pauli matrices. 
Correlation functions are obtained by repeated numerical 
integration of the classical lattice equations of motion 
and Monte Carlo sampling of initial conditions [6]. The 
initial time derivatives A° (i — 0, x) are set to zero for all 
shown results, which implements the Gauss constraint at 
all times. This simplifies the numerical implementation, 
however, we checked that including initial electric fields 
does not change much the observed growth rates. Shown 
results are from computations on = 128^ lattices. 

Using the initial condition (jj]), supplemented by a 
small noise for all gauge field components, we indeed 
find exponentially growing behavior corresponding to a 
nonequilibrium instability. In Fig. [5] we show the abso- 
lute value squared of A^{t,p) as a function of time in 
units of the fourth root of the energy density e. Dif- 
ferent lines correspond to different spatial momenta and 
one observes that modes with smallest momentum have 
the biggest initial growth rate, as expected from the lin- 
ear analysis of Sec. [TTl In contrast to the linear analysis. 
Fig. [5] exhibits that higher momentum modes can also 
have significant growth at somewhat later times. 

Concentrating first on earlier times, we obtain an av- 
erage growth rate by fitting an exponential function to 
|j4^(t,P2)| over times large compared to the oscillation 
frequency. The results of such a fit are shown in Fig. |31 
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FIG. 3. Growth rates for |j4!j(t,pz)j from simulations with 
different magnetic field strengths measured (a) in lattice units 
and (b) in units of e^''^. The solid lines represent the function 
(|30p for the various magnetic field strengths. The presence of 
a subdominant instability band is shown in the appendix, in 
FigHni 



The upper panel shows the growth rate as a function 
of pz for different initial magnetic field strengths. Both 
the growth rate and the momentum are plotted in lat- 
tice units. We find that the primary growth rates can be 
described very well by 

7^"""^^ (p.) ^ V0A2gBit^0)-pl (30) 

for momenta < 0A2gB{t = 0). This function is dis- 
played as a curve in Fig. [3] (a) along with the fit values. 
The numerical factor appearing in front of gB(t — 0) in 
(1501) is equal within errors to the one in (|27p found from 
the linear analysis. 

In the lower panel of Fig. [3] (b) we plot the same data 
as in (a), however, this time the growth rates as well as 
momenta are measured in units of {g'^eY^'^. One observes 
a collapse of the results from simulations with different 
initial field strength onto a single curve to very good ac- 
curacy. The reason for this trivial scaling can be under- 
stood from (|30l) for the growth rate and the fact that 
the energy density is proportional to B^. We note that 
parametric resonance leads to a subdominant instability 
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FIG. 4. Time evolution of the individual gauge field components in units of appropriate powers of the energy density. Each 
panel shows the time evolution of three different Fourier coefficients whose momenta are parallel to the z-axis. 



band at higher longitudinal momenta, which is given in 
the appendix. 

Fig. [2] shows that there are significant deviations from 
the linear analysis for higher momentum modes already 
at rather early times, before the growth of the lowest 
momentum mode saturates. Though these higher mo- 
mentum modes have initially comparably small primary 
growth rates, there is a substantial speed-up because of 
non-linearities at later times. As a consequence, they al- 
most catch up with the initially fastest growing mode. 
The observed non-linear behavior is the consequence of 
the self- interactions of the gauge fields [6] . The secondary 
growth rate is to a good approximation three times the 
growth rate of the fastest primary growth, which can be 
described from resummed loop expansions based on the 
two-particle irreducible effective action following the lines 
of Ref . [isj l , where similar phenomena for nonequilibrium 
instabilities in scalar field theories were studied. 

Even though the initial conditions (|3]) involve only 
few gauge field components, the non-linear dynamics 
very efficiently includes the other components as well. 
Similar to Fig. [21 we show in Fig. |4] the time evolu- 
tion for all gauge field components. The three lines in 
each graph correspond to the different momenta Pz = 0, 
Pz = 0.48(52e)i/4 and pz = f.2(g2e)i/4^ respectively. 
The zero momentum mode oscillations of AZ and A^,, 
expected from the linear analysis of Sec. are mani- 
fest as well as the primary growth in these components 
at non-zero momenta and in A^ and The simula- 
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FIG. 5. Time evolution of the diagonal entries of the spatial 
components of the stress-energy tensor, i.e. transverse ver- 
sus longitudinal 'pressure', for single runs with a macroscopic 
initial B-field configuration defined in Q. 



tions reveal growth for all other components such that 
most of the details about the initial conditions are lost 
rather quickly. All growth saturates at approximately 
t ~ 40 ((;^e)^/'*, which coincides with the time when the 
zero modes cease to oscillate. Before this time these os- 
cillations are almost identical to the unperturbed case 



^ At early times the magnetic field is approximately given by 



B\{t) ~ A^{t, 0)A^{t, 0), because the contribution of the deriva- 



tive terms is small. 
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described by ([13]). 

The initial configuration ([1} pointing in the longitu- 
dinal direction is highly anisotropic and an important 
question in this context is the characteristic time for 
isotropization. For applications to hydrodynamic de- 
scriptions of heavy ion collisions the isotropization time 
of the stress-energy tensor is of particular relevance. The 
diagonal entries of the spatial components of the stress- 
energy tensor, or transverse and longitudinal 'pressure', 
are shown as a function of time in Fig. [5j One observes 
that for the macroscopic initial i?-ficld configuration de- 
scribed by (|H) the stress-energy tensor shows oscillatory 
behavior. We find significant damping of these oscil- 
lations and approximate isotropization about the time 
when the exponential growth with characteristic rate 
\/gB stops. As we will see in Sec. IIVI this result can 
change dramatically if fluctuations are included in the 
initial conditions. 
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FIG. 6. Growth rates for the fluctuations of gauge potentials, 
{\A^{t,pz)\'^)^^'^ etc., for the ensemble of coherent fields as 
described in the text. 



IV. INCLUDING INITIAL FLUCTUATIONS 

A. Ensemble of coherent fields 

So far, we considered initial configurations with a 
sharply defined macroscopic field amplitude B described 
by ([¥]). An important step towards more realistic sce- 
narios concerns the inclusion of fluctuations in the initial 
conditions. Here we choose non-zero initial values of the 
homogeneous A° fields randomly from a Gaussian distri- 
bution with zero mean and finite width. Then we perform 
an ensemble average over the initial conditions. Such a 
nonequilibrium ensemble is meant to describe a single 
system in the same sense as, e.g., equilibrium thermody- 
namics may be applied to a macroscopic system. More 
precisely, we consider at initial time the homogeneous 
field configurations 



gA^^si, gA.y = s2, 
where the real random numbers si and S2 fulfill 



(Sl> = (S2> =0, 



(31) 



(32) 



where (. . .) denotes the ensemble average with given 
width A. Again, all other Af are taken to vanish initially 
up to a small amplitude noise seeding the instabilities. 

It turns out that growth rates in units of the average 
energy density (e) are fairly independent of the value of 
the width and the lattice size for the considered range 
of parameters. Fig. [S] shows results for growth rates us- 
ing the range of width Aa^ — 0.15 — 0.25 ensuring that 
the instability dynamics can be well resolved on the lat- 
tices we are using. We employ 64'^ or 128'^ lattices for 
an ensemble built out of 1600 runs in total. The ver- 
tical lines indicate the statistical errors. The ensemble 
averages show a rather similar picture as for the initial 
condition (j?]) without fluctuations. For the latter case 
we saw in Sec. ITT] that SAf, - SAl and SAl + SAl exhibit 
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FIG. 7. The same as in Fig. (5] however, including initial 
fiuctuations described by the ensemble of coherent fields (|31[) . 
The vertical bars reflect the size of typical fluctuations on the 
smoothed curve we plot. The thin dashed line (1/3) corre- 
sponds to an isotropic Tij. 



primary growth rates as displayed in Fig. JH] Also for the 
ensemble average the initially unstable modes turn out 
to be A^, Ay, Ay and A^. However, two different growth 

rates are observed from Fig. [5] Here (|^^(i,pz)p)^/^ and 
{\Ay{t,pz)\'^)^^^ exhibit a somewhat bigger growth rate, 
while (|A3(i,p^)|2)i/2 and (|A2(i,p^)|2)i/2 have a smaller 
rate than for the simple initial condition Q. 

The inclusion of initial fluctuations has a significant 
impact on the behavior of the stress-energy tensor. In 
Fig. [7] we plot the diagonal entries of the spatial com- 
ponents of the stress-energy tensor as a function of 
time. Even though the initial configuration is highly 
anisotropic, one observes that transverse and longitudi- 
nal pressure very quickly approach the same value. In 
contrast to what has been observed in Sec. IIIIl where ap- 
proximate isotropy of pressure occured at the end of the 
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FIG. 8. Time evolution of a low momentum mode of the FIG. 9. The growth rate of a function of the momentum for 

gauge potential Ay{t,pz — 0) for different sizes of correlated various domain sizes. 

domains. 



period of exponential growth of fluctuations, it now hap- 
pens on a much shorter time scale. This rapid isotropiza- 
tion is a dephasing phenomenon. For a sharply defined 
initial field amplitude B, which is realized with the initial 
condition the longitudinal and transverse pressure 
components oscillate with frequency ^ y/'gB according 
to Fig. [S] The time average of the pressure over several 
oscillations was already isotropic in that case, but the co- 
herent oscillations of the homogeneous A fields showed up 
in this quantity. In contrast, sampling many runs with 
different phases and initial field amplitudes leads very 
quickly to a practically constant value. Rapid isotropiza- 
tion is an important ingredient for the application of hy- 
drodynamic descriptions for heavy-ion collisions. How- 
ever, in this case longitudinal expansion and increas- 
ing diluteness may make the system again become more 
anisotropic. We also emphasize that the observed de- 
phasing phenomenon depends on the fact that the stress- 
energy tensor is dominated by the homogeneous modes 
for the considered initial conditions. For the spatially in- 
homogeneous configurations in the context of color flux 
tubes, discussed in the following, this will iro longer be 
the case. 



B. Modeling color flux tubes 

In this section, we also take into account that the fields 
may be correlated only over a limited transverse size 
while being homogeneous in longitudinal direction. To 
describe this, we divide the transverse plane into square- 
shaped domains of equal size. Each patch is filled with 
a coherent color magnetic field as in (|31|) . however, in 
each domain we use a different, randomly chosen color 
direction. The distribution of the random yl-fields is a 
Gaussian with a given width A. This corresponds to an 
inhomogeneous configuration, which is built from corre- 
lated domains of a characteristic transverse size. If the 
transverse size is taken to be Qs then this should model 



characteristic aspects of the dynamics of color flux tubes. 
The typical scale of the magnetic field inside a domain is 

^/-^domain ^ A . 

We consider domains ranging in size from 0.4 to 6 in 
units of the characteristic spatial extent 1/A. As be- 
fore, also a small amplitude noise is added to all the 
fields to seed the instabilities. Similar to the dynamics 
for the homogeneous initial conditions described above, 
one observes also for the inhomogeneous configurations 
approximately exponential growth of the spatial momen- 
tum modes of the gauge potentials, Af{t,p). However, 
one observes primary growth for all color and spatial 
components with approximately the same rate. In Fig. [S] 
the time evolution for a low momentum mod£0 of the 
gauge potential, A\{t,pz — 0), is shown as a function of 
time in units of the averaged energy density e for systems 
with different initial domain sizes. The data is obtained 
from an average over only eight runs, but since each run 
contains a significant number of domains with indepen- 
dent coherent fields it turns out that the growth rates are 
fairly stable. 

In Fig. [S] the growth rates are given as a function of 
momenta for different domain sizes. Comparing the re- 
sults for the homogeneous and the inhomogeneous initial 
configurations, one observes a similar momentum depen- 
dence, but the latter show smaller growth rates by ap- 
proximately a factor of six to eight. It is interesting 
to note that one still might understand this difference 
in terms of the characteristic size of the space-averaged 
color-magnetic fields. This average is not i^domain be- 
cause different domains have independent fields, which 
can average out. Instead, the characteristic size may be 
associated to the average of the absolute value of the 
magnetic field. In the case of the 'flux tube' initial con- 
dition the space-average of the magnetic held is typically 
a factor of hundred smaller than in the homogeneous case 
in units of the energy density. If the growth rates scale 



^ We use the first non-zero momentum mode on a 128^ lattice. 
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with the square root of the magnetic field, as is char- 
acteristic for a Nielsen- Olesen instability, one expects a 
factor of ten smaller growth rates. However, this is only 
a lower estimate, since somewhat larger growth rates can 
occur if the spatial variations of the magnetic field are 
such that there are regions with larger magnetic fields 
which cancel in the average. 

Remarkably, the growth rates measured in units of the 
energy density of the system are similar to the growth 
rates observed in Ref. j6l|. In the latter study inhomo- 
geneous initial field configurations were randomly gen- 
erated from a Gaussian distribution. The oblate initial 
distribution for the spatial Fourier modes = 0, p)| 

was characterized by a transverse width At and a lon- 
gitudinal width Al, where A^ ^ A^. This corresponds 
to an extreme anisotropy with a (5(p^)-like distribution 
in longitudinal direction. As a consequence, the approx- 
imately homogeneous field configurations in longitudinal 
direction employed in Ref. ^ resemble quite closely the 
correlated domains employed here if the transverse width 
At is associated to the inverse domain size of our con- 
figurations. 



V. CONCLUSIONS 

In this paper we have studied out-of-equilibrium dy- 
namics in SU (2) pure gauge field theory. We presented a 
detailed analytic calculation of the time evolution in the 
linear regime starting from coherent field configurations. 
Generalizing the well-known Nielsen-Olesen instabilities 
for constant color-magnetic fields, in a first step we took 
into account temporal modulations. This leads to a re- 
markable coexistence of the original Nielsen-Olesen insta- 
bility for time-averages together with a subdominant in- 
stability band because of the phenomenon of parametric 
resonance. The comparison to classical-statistical lattice 
gauge theory simulations showed remarkable agreement 
between analytics and numerics in the linear regime, 
while the lattice results show important nonlinear phe- 
nomena at later times such as secondary growth rates 
which are multiples of primary growth rates. 

This analysis provided the building blocks for the un- 
derstanding of the dynamics of more realistic initial con- 
ditions including fiuctuations. For the ensemble of co- 
herent fields, which have zero mean color-magnetic field, 
again robust growth of the Nielsen-Olesen type is ob- 
served. Remarkably, we find that for this ensemble 
isotropization of the stress-energy tensor happens on a 
much shorter time scale than the characteristic inverse 
growth rate of the instability. This rapid isotropiza- 
tion is a dephasing phenomenon. Rapid isotropization 
is an important ingredient for the application of hydro- 
dynamic descriptions for heavy-ion collisions. However, 
in this case longitudinal expansion and increasing di- 
luteness Q can make the system again become more 
anisotropic. In particular, the observed dephasing phe- 
nomenon depends on the fact that the stress-energy ten- 



sor is dominated by the homogeneous modes for the con- 
sidered initial conditions. For spatially inhomogeneous 
configurations this is typically not the case. 

To see this and in order to make the link to earlier 
work on plasma instabilities using inhomogeneous initial 
field configurations, we finally took into account that the 
fields should be correlated only over a limited transverse 
size. If the transverse size is taken to be Qs or ^/gB then 
this should model characteristic aspects of the dynamics 
of color fiux tubes. In this case we observed exponential 
growth for all components of the inhomogeneous gauge 
field configuration. The nature of the underlying insta- 
bility still shows properties reminiscent of the Nielsen- 
Olesen-type. In particular, the results indicate that the 
characteristic dependence of the primary growth rate as 
a function of momentum still reaches its maximum at 
low momenta. In contrast, Weibel instabilities are ex- 
pected to show a vanishing growth rate at zero momen- 
tum , which has also been emphasized in Ref. ilOj . In 
view of these results it is interesting to observe that pre- 
vious findings about plasma instabilities from classical- 
statistical lattice gauge theory |^, Q share important 
aspects of a Nielsen-Olesen instability. 

In this work, we have not discussed the extraction of 
distribution functions and their possible interplay with 
time-dependent condensates after isotropization as sug- 
gested recently in Ref. |23||. Distribution functions may, 
for instance, be derived from equal-time correlation func- 
tions in Coulomb gauge using the classical-statistical sim- 
ulations employed in this work. This is deferred to a 
separate publication 2^ . 
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APPENDIX 

A. Linear regime including parametric resonance 

In Sec. |lT]we discussed analytic solutions in the linear 
regime, which concentrated on the dominant exponen- 
tially growing modes at low momenta. However, the 
analysis neglects the phenomenon of parametric reso- 
nance in the presence of a periodically evolving back- 
ground field. In this appendix we show that this phe- 
nomenon leads to a subdominant instability band at 
higher longitudinal momenta than the Nielsen-Olesen in- 
stability. Here we focus on the evolution of 6A-{t,p) 
described by ([2^ of Sec. |lT]as this shows the dominant 
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primary instability. Plugging in the solution for the evo- 
lution of the background field p^ . the evolution equa- 
tion for vanishing transverse momentum reads 



d 



- gB CT? 



/gS t, - 



5A^{t,p) = 0,(33) 



where in this appendix we always denote the longitudinal 
momentum by p to ease the notation. This closely resem- 
bles the Jacobian form of the Lame equation. The crucial 
difference here is the negative sign in front of the oscil- 
lating term, which gives rise to the Nielsen-Olesen type 
instability discussed in Sec. [TTl Our strategy for solving 
the above evolution equation consists of approximating 

cn^(a;, m) = 1 - cn^(.T - K{m),m) + 0{m^) (34) 

for TO = 1/2 in ([55)) in order to reduce it to the well 
known Lame equation, where again K(m) denotes the 
complete elliptic integral of the first kind [1^ . While this 
approximation captures the important features of (j33p , it 
has the shortcoming that the average field strength gB 
according to Eq. (P?]) is overestimated by a factor of 



r(3/4)2 



1 



r(5/4)r(i/4) 



0.457 . (35) 



Consequently, we expect the approximation to yield 
slightly enhanced growth rates that extent to somewhat 
higher momenta within about ten percent accuracy of 
the full numerical solution at early times, where the lin- 
ear analysis is expected to hold. In order to obtain the 
solutions for SA-{t,p) for the above approximation, we 
recast the evolution equation to the Weierstrass form by 
use of the identity [11 |2l| 



cn 



1 

"'2 



= -2p 



(36) 



where p{x) is the Weierstrass elliptic function with roots 
ei — 1/2, 62 = and 63 = —1/2. The evolution equation 
in Weierstrass form then reads 



P 



gB - 2gB p ( ^5^* - r))] M_ {t,p) = , 

(37) 



where r = (1 — i)K{l/2). By introducing the dimen- 
sionless time variable 9 — \fgB t and expressing the time 
independent contribution in (jSTj in terms of the Weier- 
strass elliptic function according to 



El 
gB 



(38) 



the fundamental solutions Up{6) and Up{9) to ([57)1 read 



„T(e-T)c(z) 



a{e-T± z) 
a{9-T) 



(39) 



Here the superscript 1/2 labels the respective fundamen- 
tal solution. The functions (t{x) and Ci^) represent the 
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numerical 
rescaled analytic 

Yresc=T(Pz(Bapp/B)^")(B,pp/B)"2 




FIG. 10. Growth rates and oscillation frequencies in the linear 
regime. 



corresponding Weierstrass functions and the momentum 
dependence of the solutions is encoded in z defined by 



In order to obtain the growth rate ^(p) and the oscilla- 
tion frequency ui{p) we proceed along the lines of Sec. [ID 
performing a Floquet analysis, i.e. 



_,^cy^(p)c/i/2(0) (40) 



U^^^9 + 2K{l/2)) = e 

with the Floquet function Cy^(p) as in SecHTj Using the 
quasi periodicity of the Weierstrass cr-function [l^, [2lj 

a{x + 2K{l/2)) = -a{x) exp[2(a; + K{l/2))(l{K{l/2))] 

(41) 

we find 

c^j\p) - ±2* [i^(i/2)c(z) - z amim ■ m 

The growth rate ^(p) and oscillation frequency w(p) are 
related to the real and imaginary parts of the floquet 



index b'v 



l{p) 
u}{p) 



gB 
'gB 



\MC-{p)]\ 

2K{l/2) 

\MC-{p)]\ 
2K{l/2) 



(43) 
(44) 



^ We note that in this way the oscillation is only defined up to 
constants of 27r. 



11 



In order to evaluate the expressions in (|43)) and (1441) one 
needs to evaluate the mapping ([55]) defining z. Follow- 
ing the lines of [2l[ we find that there are four different 
regimes described by 



with 



2 = /3 



Z = ;9 + iif(l/2) 



i/3 



for — - 



gB ^2 



(45) 



for 



for 



for 



^^3 
55 2 



(48) 



where /3 S [0, if (1/2)] for all regimes. The different 
regimes correspond to two stable and two unstable bands. 
The lowest band corresponds to the Nielsen- Olesen insta- 
bility, whereas the third band corresponds to the para- 
metric resonance instability. Outside these momentum 
regions all modes are stable within the approximation 
using ((M)) . 

The corresponding growth rates and oscillation fre- 
quencies for all bands are on display in Fig. [TU] The 
modes being subject to the Nielsen- Olesen instability 
show no oscillatory behavior according to this approx- 
imation. In contrast, modes which are amplified by 
the parametric resonance instability follow the oscilla- 
tion of the macroscopic field A(t). The stable band in 
between interpolates between the two oscillation frequen- 
cies, while at higher momenta ^ gB the free field limit 
'^(p) — V is approached. Concerning the growth rates 
one observes that the Nielsen- Olesen instability exhibits 
the dominant growth rate when compared to paramet- 
ric resonance. When comparing the analytic results for 
the growth rate to lattice data, we find that the ana- 
lytic prediction is somewhat larger and extends to higher 
momenta. As mentioned above, this is a consequence 
of the approximation p4p . However, one can accurately 
describe the lattice data by a trivial rescaling of the mag- 
netic field amplitude B in the analytical result by a factor 
of \/c/{l — c) as suggested by ((35|) to reproduce the cor- 
rect average amplitude B. 



B. Parametric resonance band 

The above discussion provided implicit expressions for 
the growth rate as a function of momentum. To obtain 
explicit expressions one has to evaluate the variable z 
introduced in ([38]) as a function of the momentum p. 
According to (|45l) -(|48 l) this has to be done separately 
for each band. In Sec. |TI] we have already provided an 
explicit expression for the growth rate of the Nielsen- 
Olesen instability by an alternative approach. The same 
is desireable for the parametric resonance instability. To 
obtain such an expression we evaluate z as a function of 
p, where for the parametric resonance band (1471) one has 



z ^ (3 + iK{l/2) 



(49) 




(50) 



The last expression was obtained by expressing the 
Weierstrass function p(z) in (p8)) in terms of Jacobi el- 
liptic functions according to Though ([50]) along 
with (j43p already provides an explicit expression for the 
growth rate, it proves insightful to apply further approx- 
imations. In a first step we expand the Weierstrass zeta 
function as 



k=l ^ 



. = C(-) = ^ + 0(.^) 



sm 



/ knx\ 

\ UJ J 

(51) 
(52) 



where q = e ^ and oj ~ K{l/2). Without further ap- 
proximations the growth rate 7(73) is then given by 



lip) 



Re 



TT / TTZ \ 

t;- cot — 
2uj \2ujJ 



2ti 



UJ ^ 1 — Q 

fc=l ^ 



/ k'Kz\ 

2^ "° \— ) 



(53) 



The real part can then be expanded in a q-Series by 
seperately expanding the cotangent and the series of sine 
functions. As g ~ 0.043 is a reasonably small expansion 
parameter we will keep the leading terms only. For the 
cotangent the real part is given by 



Re 



TT 

cot I —p + i — 
\2uj^ 2 



sin(5/3) 



cosh(7r) - cos(^/3) 
= 2gsin(^/3) +0(9^) , (54) 

and similarly for the series of sine functions one finds 



4Re 



00 

.fc=i 



^UJ 
Jik 



Bin + ikn)j 

in 

2(7sin(^/3) +0{q^) 



n2k 



00 u 
k=l 



(55) 



As these expressions involve trigonometric functions of 
f3(p) a simple expression for the momentum dependent 
growth rate 7(p) in the parametric resonance regime can 
be obtained by expanding the inverse of the Jacobi co- 
sine in (1501) in terms of inverse trigonometric functions 
according to 



cn 



'^{x, m) 



2uj 



■ cos 



-^{x)+0{m) . 



(56) 
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By use of double angle formulas we obtain as the final 
result 



for the parametric resonance band where 1 < j (gB) < 
3/2. The maximal growth rate 70 is realized at the mo- 
mentum p ~ ^Jh]\\fgB . Numerically one finds 70 ~ 
0.146-^^5 which is significantly smaller than the growth 
rates associated with the Nielsen-Olesen instability. 
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